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Abstract 

We propose a non-perturbative procedure to fix generic covariant gauges 
on the lattice. Varying the gauge parameter, this gauge fixing provides a 
concrete method to check numerically the gauge dependence of correlators 
measured on the lattice. The new algorithm turns out to converge with 
a good efficiency. As a preliminary physical result, we find a sensitive 
dependence of the gluon propagator on the gauge parameter. 



1 Introduction 



Lattice gauge fixing is unavoidable to compute non perturbatively the propaga- 
tors of tlie fundamental fields appearing in the QCD Lagrangian [^. It is also 
necessary in non-perturbative renormalization schemes p|, ^ which use gauge de- 
pendent matrix elements to renormalize composite operators, and it becomes a 
fundamental technical ingredient in the so called non gauge invariant quantiza- 
tions of chiral theories @]. 

Up to now, the Landau gauge is the only covariant gauge for which an efficient 
numerical algorithm is available. 

A few years ago, an unconventional method to implement a generic covariant 
gauge on the lattice was proposed 0. In principle this procedure takes into 
account the contributions from each gauge orbit with appropriate weight. With 
this method a gauge dependence of the gluon propagator, in particular at zero 
momentum, had been found On the other hand the implementation of this 
method for physical lattices can be numerically demanding. 

Here we follow a more conservative procedure to fix a generic covariant gauge 
on the lattice that stems from ref. [|7|. It is based, as in the Landau case, on the 
minimization of a functional Ha[G\ chosen in such a way that its absolute minima 
correspond to a gauge transformation G satisfying the appropriate gauge condi- 
tion (see section In the continuum it is easy to show that this procedure 
is equivalent to the Faddeev-Popov quantization for covariant gauges, i.e. the 
two procedures lead to the same matrix elements for a generic gauge dependent 
operator. As usual the Faddeev-Popov factor can be written as a Gaussian inte- 
gral of local Grassman variables, the resulting effective action is invariant under 
the BRST transformations and the correlation functions of the operators satisfy 
the appropriate Slavnov- Taylor identities. On the compact lattice the presence 
of Gribov copies is an obstruction to define partition functions with exact BRST 
invariance as it has been shown by H. Neuberger a few years ago [Q. The Neu- 
berger's argument concerns a fundamental theoretical point whose relevance for 
the measurement of gauge dependent operators on the lattice must be clarified. 
Nevertheless, we believe that new numerical gauge fixing algorithms represent a 
concrete tool to improve the understanding of the lattice gauge properties. 

In the following we propose a new numerical gauge fixing procedure which 
generalizes the usual Landau gauge fixing to generic covariant gauge. The nu- 
merical implementation of our procedure on the lattice uses a straight-forward 
generalization of the standard Landau algorithm, i.e. a steepest descent iterative 
algorithm minimizes a discretization of the new functional Ha[G\. In order to 
have an efficient and feasible gauge fixing procedure, the discretized functional 
Hjj[G] must be chosen carefully. In fact a naive discretization of this functional, 
due to its complicated structure, would lead to an algorithm which either does 
not converge or takes too much computer time |]TU| . 



The main purpose of this paper is to show that a simple discretization Hjj[G] 
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exists and leads to a minimization algorithm which converges with a good effi- 
ciency. The application of this algorithm to the non-perturbative evaluation of 
gauge dependent Green functions like the quark and gluon propagator will clarify 
the gauge dependence of the fitted parameters. 

In order to show the feasibility of our method, we have computed the gluon 
propagator at different values of the gauge parameter at small volumes, finding 
a sensitive gauge dependence. 

The plan of the paper is as follows. In section |^ we review the method 
proposed to fix a generic covariant gauge on the lattice. In section || we describe 
a very simple and efficient discretization of the gauge fixing functional. In section 
^ we give the details of our numerical simulations and report our main numerical 
results. 



2 Covariant Gauge Fixing 

In this section we set the notation and briefly formulate the covariant gauge fixing 
method in the continuum. In the following we will neglect the problem of Gribov 
copies, assuming that a gauge section in the space of gauge fields intersects all 
gauge orbits once and only once. In the Landau gauge the expectation value of 
a gauge dependent operator is given by 

(O) = J M^(5ry5^C»e-^(^)-^«''-'('''*''^)5(a^A^) , (1) 

and the gauge fixing condition can be enforced non-perturbatively by minimizing 



Gribov's functional JlT 

FAG] ^ \\A^\\' = / Tr (Aj/lJ) rf^x . (2) 

The Landau gauge is readily extended to a general covariant gauge-fixing condi- 
tion of the form 

d.A^^ix) = Aix) , (3) 

where A{x) = A"(x)^ belongs to the Lie algebra of the group and Tr(T"r'') = 
26"-^. Since gauge-invariant quantities are not sensitive to changes of gauge con- 
dition, it is possible to average over A{x) with a Gaussian weight 

{0) = J SAe~^ / '^'-^'•[^'] J SA^6riSf] O e-^(^)-^s''-*('''*''^)5(5^A^ - A) , (4) 

obtaining the standard formula 



This formula is adopted as a definition for the expectation value of gauge de- 
pendent operators. The invariance under BRST transformations of the effective 
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action in eq. leads to the Slavnov- Taylor identities among different correlation 
functions. In particular the longitudinal part of the gluon propagator has to be 
equal to the free one, i.e. 

i:i{d^A^{x)d,A,{y)) oc ad{x - y) , (6) 
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where N is the number of colors. Following the usual technique, the gauge-fixing 
condition (§) is obtained non perturbatively by minimizing the new functional [0 

Ha[G] ^ J d'xTi [{d,A^ - A)(a,A^ - A)] , (7) 

which obviously reaches its absolute minima {Ha[G] = 0) when eq. @ is satis- 
fied. Therefore in this case the Gribov copies of the equation are associated 
with different absolute minima of eq. (^. All the stationary points of Ha[G] 
correspond to the following gauge condition 

^.a,(a^Aj-A)=0, (8) 

where Di, is the covariant derivative. Hence, due to the complexity of the new 
functional, its minimization could find "spurious" solutions which correspond 
to zero modes of the operator D^di, and do not satisfy the gauge condition in 
eq. (^. Of course the numerical minimization of the discretized version of eq. (0) 
can reach relative minima (spurious solutions) with Ha[G] which are not 
distinguishable from the absolute minima. Hence this could simulate the effect 
of an enlarged set of numerical Gribov copies. Within the limitations of this 
preliminary study we do not find any practical difference between the use of the 
new functional with respect to the standard Landau one (see results given in 
Fig. 1^ below). 

On the lattice, the expectation value of a gauge dependent operator O in a 
generic covariant gauge is 

(O) = ly rfAe-^E.^KA^l J d[/0(t/^'^)e-^^(^) , (9) 

which is the straight-forward discretization of Eq. (^), where A is dimensionless 
on the lattice. S{U) is the Wilson lattice gauge invariant action and Ga is the 
gauge transformation that minimizes the discretized version of the functional (0). 
On the lattice, the correct adjustment to the measure is included in eq. (||) by 
evaluating the operator over the gauge rotated links. Therefore it is not necessary 
to introduce ghost fields but it is mandatory to fulfill the gauge fixing condition 
numerically in order to get 
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3 The Driven Discretization 



In the Wilson discretization of gauge theories, the fundamental fields are the 
links which act as parallel transporters of the theory. Hence the lattice fields 

are derived quantities which tend to the continuum gluon field as the lat- 
tice spacing vanishes. As a consequence on the lattice it is possible to choose 
different definitions of formally equal up 0(a) terms. In quantum field the- 
ory this ambiguity is well understood because any pair of operators, differing 
from each other by irrelevant terms, will tend to the same continuum operator. 
This feature, checked in perturbation theory, has been verified numerically at the 
non-perturbative level in ref. ||12|, where it has been shown that different defini- 
tions of the gluon field give rise to Green's functions proportional to each other, 
guaranteeing the uniqueness of the continuum gluon field. 

The freedom to choose the lattice definition of can be used to build dis- 
cretized functionals which lead to efficient gauge-fixing algorithms. In the stan- 
dard Landau gauge fixing, for example, the discretization of the functional (0) is 
given by 

1 



Fu[G] 



VTa?g 



(10) 



where V is the 3-dimensional volume and T the time size of the lattice. This 
formula corresponds, only up to 0{a) terms, to the naive discretization of eq. (|^) 
that is obtained from the standard lattice definition of the gluon field 



A^{x) 



2iag 



Traceless 



On the other hand -FVfC] has the important property that it depends only linearly 
on G{x), when the iterative algorithm visits the lattice point x. This feature 
would be spoiled if F(/[G'] were defined assuming literally the naive discretization 
([TT]) in eq. (0). 

In order to study the convergence of the algorithm, two quantities are usually 
monitored as a function of the number of iteration steps: ^(/[G] itself and 



VT 



Y,Tr[ApA 



(12) 



where 



Ap(x) 



Xf{x) - X^x) 



Traceless 



OC 



S Fu[G] 
6e 



(13) 



being G = e"^" 2 and 



Xf{x) 



E 



U^,ix] 



+ Ul{x 



(14) 
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Ap{x) is proportional to the first derivative of F;7[G] and readies zero as the func- 
tional is extremized. Therefore the quality of the gauge fixing is determined by 
the parameter 6f which corresponds to the continuum quantity / (i^xTr(9^74^)^. 

A naive discretization of Ha[G] will generate a quadratic dependence on G, 
which could prevent the convergence of the algorithm. This obstacle has been 
overcome by taking advantage of the freedom to choose the gluon field definition. 
In fact, as in the Landau case, it has been possible to find a discretization of 
Ha[G] ("driven discretization") that depends linearly on G{x) and corresponds, 
up to 0(a) terms, to the continuum limit, eq. (0). This aim can be reached 
by choosing each different term of Hu[G] in order to guarantee the local linear 
dependence on G{x), instead of being the algebraic consequence of a particular 

definition. We propose the following compact form of Hu[G]: 

Hu[G] = y^,TrY.J''{^)J''\^) , (15) 

where 

J{x) = N{x) - igA{x) , 

N{x) = -8I + Y.(ul{x-u) + U,ix)) . (16) 



It is easy to see that locally Hu[G] transforms linearly in G{x) and its continuum 
limit is the functional (^. Hu[G] is positive semi definite and, unlike the Landau 
case, it is not invariant under global gauge transformations. 
The functional Hu[G] can be minimized using the same numerical technique 
adopted in the Landau case. In order to study the convergence of the algorithm, 
two quantities can be monitored as a function of the number of iteration steps: 
Hu[G] itself and 

Oh = ^T^TtIAhAU , (17) 



where 

Ah{x) = \Xh{x)-XIj{x) 



SHu[G] 

oc r 

Traceless 



and 

Xnix) = ^{u^{x)J{x + ^) + Ul{x - ^i)J{x - ^) 

- SJ{x) -721 + igN{x)k{x) . (19) 

Ah is the driven discretization of the eq. (|^); it is proportional to the first deriva- 
tive of Hjj[G] and, analogously to the continuum, it is invariant under the trans- 
formations A(x) — A(x) + C, where C is a constant matrix belonging to the 
5'[/(3) algebra. During the minimization process Oh decreases to zero and Hjj[G] 
becomes constant. The quality of the convergence is measured by the final value 
of Oh. 
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4 Numerical Simulations and Results 



We have generated 50 SU{3) thermalized link configurations using the Wilson 
action with periodic boundary conditions at (3 = 6.0 for 8^ and 8^ x 16 volumes. 
Following the prescription contained in eq. @, for each Monte Carlo configuration 
and for each lattice site we have extracted a matrix A(x) according to a Gaussian 
distribution at a fixed a value. The gauge-fixing code implements an iterative 
overrelaxed minimization algorithm for -F(/[(j'] and Hu[G]. We have monitored 
the quantities Fu and Op for the standard Landau gauge-fixing algorithm and 
Hu[G] and 9h for the new one after every lattice sweep. We have found that the 
value of the overrelaxing parameter u adopted for the Landau algorithm u = 1.72 
is a good choice also for the new gauge fixing. 

In Fig. [1| we report the values of 9h for different a values as a function of 
the gauge fixing sweeps for a typical thermalized configuration at /? = 6.0 and 
V ■ T = 8'^. Each sweep of the new algorithm takes ~ 15% more computer 
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Figure 1: Typical behaviour of 9h vs gauge fixing sweeps for different choices of a, for 
/? = 6.0 and F • r = 8l 



time with respect to the Landau case. For a = 0, which in this new procedure 
corresponds to the Landau gauge and is obtained by taking A = 0, the number of 
sweeps necessary to fix the gauge is ~ 1.2 times that of the standard algorithm. 
In the case of a 7^ 0, the number of sweeps to minimize Hu[G] with a given 
precision increases when a decreases. 
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Once the configuration have been rotated in a given gauge, we have tested 
our procedure computing the following two point correlation functions 



{AoAo){t) = -^^Tr(Ao(x,t)Ao(y,0)) , (20) 

{A,A,)it) ^ -i^5:5:Tr(A,(x,t)A(y,0)), (21) 

{dAdAm ^ l-j:Y.M9,M^^t)d.M^,0)) , (22) 

where and z/ run from 1 to 4, z from 1 to 3 and the trace is over the color 
indices. In eqs. (pOD-(p2D the definition (|TT|) is adopted for the gluon field on 



the lattice and indicates the usual backward derivative. The correlators in 



eqs. (^) and ( plD are relevant to the investigation of the QCD gluon sector, 
while we use the correlation in eq. (0) to check the Slavnov- Taylor identity for 
the longitudinal component of the gluon propagator. The statistical errors for the 
correlation functions have been estimated by the jacknife method. A thorough 



study of these operators will be presented in a forthcoming paper ||T3 . 

In Fig. H the Green functions computed with the standard Landau gauge 
fixing {AiAi)F and {AiAi)H evaluated with the new algorithm with a = 0, are 
reported. The remarkable agreement shows that the new algorithm reproduces 
very well the results of the standard one within an overall scale factor close to 1. 
Moreover, the possible spurious solutions of the new gauge condition apparently 
do not affect this operator at least at a = 0. 

Although in this paper we did not perform a systematic study of the Gribov 
copies for this gauge fixing procedure, a rough preliminary search shows, for 
a = 0, the same pattern of copies of the Landau case. 

In Fig. 1^ the data for the correlation (^0]) are reported. For a = they show 
the well known flat behaviour forced by the Landau gauge. For a the data 
show an almost flat curve with an enhancement at t = 1 probably due to the effect 
of contact terms. In this case the flatness stems from the fact that / d^xdoAo ~ 
because the average value of A on each time slice is negligible. 

In Fig. Hour results for the correlations {AiAi)H for different a values as a 
function of the time slice t are reported for a set of 50 SU(3) configurations at 
P = Q.O andV -T = 8\ 

Although in the case of small volume and a small number of configurations, 
the gauge dependence of the gluon propagator is clearly shown. For increasing 
a values, the time dependence of the gluon propagator becomes flatter since it 
approaches a limit where the gauge fixing effects disappear. It is also interesting to 
note that the a dependence of the gluon propagator shown in Fig. ^ does not seem 
to be re-absorbed by an overall scaling factor. This plot shows the feasibility of 
this procedure to study the gauge dependence of physically interesting correlators. 
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Figure 2: Comparison between {AiAi)F (see eq. (^l|)) obtained with the old functional 
(circles) and {AiAi)H obtained with the new one Hi/[G] with A = (crosses) as 
function of time t for a set of 50 thermalized configurations at /3 = 6.0 with a volume 
V ■ T = 8^. The errors are jacknife and the data have been slightly displaced in t for 
clarity. 



We complete the numerical checks of the gauge-fixing we are proposing, show- 
ing the numerical behaviour of a particular Slavnov- Taylor identity. This identity 
is a direct consequence of our way to implement the covariant gauge fixing and it 
must be understood as a numerical check of the procedure. This check measures 
the correlator of the operator dA with itself, and if this operator has been well 
gauge-fixed, it is related to the autocorrelator of the A's and then it must be a 
delta function. In fact, for a finite volume V ■ T with periodic boundary condi- 
tions, the naive Slavnov- Taylor identity gives for the correlator of eq. (|2^) the 
following relation: 

{dA OA) (t) oc a{6t,o - ^) + 0{a) . (23) 

Eq. ( p3D shows that the shape of the correlation function ( |2^ is expected to 
be proportional to a 5-function up to the square of the df^Af^{x) renormalization 
constant, apart from volume and 0{a) effects. In order to verify eq. (PHj), in Fig. ^ 
we plot the correlation functions (p2|) versus the time t. We believe that the good 
agreement between these data and eq. (^) is an indication of the absence of 
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Figure 3: Behaviour of {AoAo)h (see eq. (|2Cl|)) obtained with the new functional Hu[G] 
for different a values as function of time t for a set of 50 thermalized configurations 
at /3 = 6.0 with a volume V ■ T = 8"^. The errors are jacknife and the data have been 
slightly displaced in t for clarity. 



possible spurious solutions induced by our method. The small signal at t = 2,8 
is likely due to the residual correlation induced by the enhancement in the point 
t = 1. 

5 Conclusions 



In this paper we have described an efficient method to fix non-perturbatively 
a generic covariant gauge on the lattice. This procedure is equivalent to the 
Faddeev-Popov quantization for covariant gauges. Therefore it can be used to 
compare the lattice renormalized correlation functions with the same quantities 
computed in continuum perturbation theory. 

We have tested the algorithm on different SU(3) lattices for different values of 
the gauge parameter a. Our numerical data show that the statistical fluctuations 
of measured quantities is comparable to the standard Landau gauge fixing algo- 
rithm. 

We have computed the correlation functions relevant for the investigations of the 
gluon propagator at a few a values and we have found a sensitive dependence on 
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Figure 4: Behaviour of (see eq. (pl|)) obtained with the new functional i?;7[G] 

for different a values as function of time t for a set of 50 thermalized configurations 
at /3 = 6.0 with a volume 1/ • T = 8^. The errors are jacknife and the data have been 
slightly displaced in t for clarity. 



the gauge parameter. 

In collaboration with the Boston University Center for Computational Science, 
we are applying this method to the study of the gluon propagator on physical 
lattices. 
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Figure 5: Behaviour of {dAdA)H (see eq. (P2[) ) obtained with the new functional -ff(7[G] 
for different a values as function of time t for a set of 50 thermalized configurations at 
P = 6.0 with a volume V -T = 8^. The errors are jacknife. 
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